<!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.0 Strict//EN" "http://www.w3.org/TR/xhtml1/DTD/xhtml1-strict.dtd"><html xmlns="http://www.w3.org/1999/xhtml"><head><title>R: A pedigree data on 282 animals deriving from two generations</title>
<meta http-equiv="Content-Type" content="text/html; charset=utf-8" />
<link rel="stylesheet" type="text/css" href="R.css" />
</head><body>

<table width="100%" summary="page for meyer"><tr><td>meyer</td><td style="text-align: right;">R Documentation</td></tr></table>

<h2>A pedigree data on 282 animals deriving from two generations</h2>

<h3>Description</h3>

<p>A data frame attributed to Meyer (1989).
</p>
<p>&ldquo;The pedigrees for each of these 282 animals derive from an additional 24 base population
(Generation 0) animals that do not have records of their own but, nevertheless, are of interest with respect to
the inference on their own additive genetic values. Furthermore, it is presumed that these original 24 base 
animals are not related to each other. Therefore, the row dimension of u is 306 (282+24).&rdquo; (Templeman \&amp; Rosa 2004)
</p>


<h3>Usage</h3>

<pre>data(meyer)</pre>


<h3>Format</h3>

<p>A data frame containing 306 records</p>


<h3>Source</h3>

<p>Meyer K (1989). Restricted maximum likelihood to estimate variance components for animal models with
several random effects using a derivative-free algorithm. Genetics, Selection, Evolution 21:317-340.
</p>
<p>Tempelman RJ, Rosa GJM. Empirical Bayes Approaches to Mixed Model Inference in Quantitative Genetics.
in Saxton AM (Ed). Genetic Analysis of Complex Traits Using SAS, chapter 7. SAS Institute Inc., Cary,
NC, USA, 2004
</p>


<h3>Examples</h3>

<pre>
## Not run: 
library(gap)
meyer &lt;- within(meyer,{
   g1 &lt;- ifelse(generation==1,1,0)
   g2 &lt;- ifelse(generation==2,1,0)
})
lm(y~-1+g1+g2,data=meyer)
library(MCMCglmm)
m &lt;-MCMCglmm(y~-1+g1+g2,random=animal~1,pedigree=meyer[,1:3],data=meyer,verbose=FALSE)
summary(m)
plot(m)   

meyer &lt;- within(meyer,{
   id &lt;- animal
   animal &lt;- ifelse(!is.na(animal),animal,0)
   dam &lt;- ifelse(!is.na(dam),dam,0)
   sire &lt;- ifelse(!is.na(sire),sire,0)
})
# library(kinship)
# A &lt;- with(meyer,kinship(animal,sire,dam))*2

A &lt;- kin.morgan(meyer)$kin.matrix*2

library(regress)
regress(y~-1+g1+g2,~A,data=meyer)
prior &lt;- list(R=list(V=1, nu=0.002), G=list(G1=list(V=1, nu=0.002)))
m2 &lt;- MCMCgrm(y~-1+g1+g2,prior,meyer,A,singular.ok=TRUE,verbose=FALSE)
summary(m2)
plot(m2)   

## End(Not run)
</pre>


</body></html>
